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We present and analyse new _R-band frames of the gravitationally lensed double quasar 
FBQ 0951+2635. These images were obtained with the 1.5m AZT-22 Telescope at 
Maidanak (Uzbekistan) in the 2001—2006 period. Previous results in the R band 
(1999—2001 period) and the new data allow us to discuss the dominant kind of mi- 
crolensing variability in FBQ 0951+2635. The time evolution of the flux ratio A/B 
, does not favour the continuous production of short-timescale (~ months) flares in 

■ the faintest quasar component B (crossing the central region of the lensing galaxy). 

Instead of a rapid variability scenario, the observations are consistent with the exis- 
tence of a long-timescale fluctuation. The flux ratio shows a bump in the 2003—2004 
period and a quasi-flat trend in more recent epochs. Apart from the global behaviour 
I of A/B, we study the intra-year variability over the first semester of 2004, which is 

reasonably well sampled. Short-timescale microlensing is not detected in that period. 
Additional data in the i band (from new i-band images taken in 2007 with the 2m 
k><( [ Liverpool Robotic Telescope at La Palma, Canary Islands) also indicate the absence 

'jj ■ of short-timescale events in 2007. 

fn ' 

Key words: gravitational lensing - galaxies: general - quasars: general - quasars: 
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1 INTRODUCTION 

FBQ 0951+26 3 5 w as discovered a decade ago by 
ISchechter et all l| 19981 ). This is a double quasar (consist- 
ing of two components A and B) at redshift Zs = 1.246, 
which is gr avitationally lensed b y an early-type galaxy at 
Zi = 0.260 (|Eigenbrod et al]|2007h . The main lensing galaxy 
robably belongs to a group of galaxies at similar redshift 



Williams et al 



20061 ). 
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Optical follow-up of FBQ 0951+2635 has been done 
in the current decade. Th is includes early ima ging and 
monitoring with the HST (|Kochanek et al. |2000|) and th e 
Nordic Optical Telescope (NOT. Ijakobsson et all 120051 ). 
as well as imaging from the Sloan Digital Sky Survey 
(SDSS. lAdelman-McCarthy et al.|[200l ). The 2.5-year mon- 
itoring campaign at the NOT focused on the R passband. 
These _R-band frame s allowed Jakobsson et al.. (2005 ) and 
iParaficz et al.l l|2006l ) to study the time delay between com- 
ponents and early extrinsic variabili ty. The time delay be - 
tween A and B is of about two weeks (| Jakobsson et al.]|2005l ). 
and there is clear evide nce for f?-band extrin s ic variations in 
the 19 99-2001 period l|Paraficz et al.ll2006l ). iParaficz et al] 
(j2006l ) showed a possible gradient of about 0.1 mmag day"'^, 
as well as a possible 50-mmag event on a timescale of several 
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Table 1. New optical observations of FBQ 0951+2635. 



Observatory 


Telescope 


Camera 


Filter 


Period 


Nights 


Exposures / night 


Maidanak (Uzbekistan) 


1.5m AZT-22 


BroCam 


R 


April 2001 


1 


2 (120 + 180 s) 










March 2002 


3 


2-4 (180/240 s) 










Apr-May 2003 


2 


2 (180 s) 










Jan-May 2004 


14 


3-5 (180/210 s) 










Dec 2004-Apr 2005 


6 


4-8 (180 s) 










Nov-Dec 2005 


3 


2 (300 s) 










Apr-May 2006 


8 


10 (180 s) 


La Palma (Canary Islands) 


2m Liverpool 


RATCam 


i 


Feb-May 2007 


52 


5 (100 s) 



months. These extrinsic fluctuations (which are not origi- 
nated in the source quasar) were attributed to microlens- 
ing by c ollapsed objects within the lensing galaxy (e.g., 
IWambsgan ss 1990, and references therein). 

However, the information obtained in the first years of 
monitoring with the NOT (1999-2001) does not permit to 
decide on the true microlensing variability, so additional 
imaging is required to address this issue. For example, i n 
Fig. 2 (middle and bottom panels) of iParaficz et al.l (|200d ). 
one sees that the microlensing gradients can account for ba- 
sically all variability without the need of introducing addi- 
tional short-term (~ months) microlensing events (see the 
distributions of points around the linear fits). Alternatively, 
in the same panels of Fig. 2 of iParaficz et all ((20061), 

sets of 

short-term microlensing events can also explain the observed 
variations. Thus, it is unclear what kind of microlensing fiuc- 
tuations occur in FBQ 0951-1-2635: short-timescale events or 
gradients lasting years (tracing a long-timescale event)?. Al- 
though both kinds of fiuctuations may be present, one rea- 
sonably expects the presence of a dominant kind accounting 
for most of the observed microlensing variability. In this pa- 
per, we try to identify the dominant microlensing variations. 

In Sect. 2 we analyse new 7?-band images taken at the 
Maidanak Observatory in the 2001—2006 period. The pre- 
vious A^Or results and the Maidanak analysis permit us to 
check the kind of microlensing fluctuations in the red part 
of the optical continuum (in the observer frame). In Sect. 3 
we study FBQ 0951-1-2635 by using a redder optical filter 
{i Sloan passband). The i-band frames correspond to a new 
monitoring campaign with the 2m Liverpool Robotic Tele- 
scope (LRT) during the first semester of 2007. To support 
our main results, we also use a few complementary frames in 
public and private archives. The conclusions are presented 
in Sect. 4. 



2 LIGHT CURVE AND FLUX RATIO IN THE 
R BAND 

FBQ 0951+2635 is part of a compact lens system, 
since the two quasar components are separated by I'.'l 
l|Schechter et al.l [19981 ) and the very faint lensing galaxy is 
0'.'2 a way from the faintest component B (jjakobsson et al] 
I2OO5I ). Fortunately, the lensing galaxy remains too faint to 
be detected with a standard R filter ( Jakobsso n ct al...2005< ). 
and this is an advantage when doing photometry. Thus, the 
system can be described by two stellar-like objects. There 
are two field stars near the lensed quasar: the bright star 
SI and the faint star S3 (see Fig. The FWHM of the 



seeing disc is measured on the SI star, which is also used 
to estimate the point spread function (PSF) of the stellar- 
like sources. This estimation allows us to obtain PSF fitting 
photometry for the double quasar. 

A single-epoch magnitude difference (ms — m^) could 
not represent the magnitude difference at the same emis- 
sion time (Amsyi), because one should t ake into account 
the 1 6-day time delay between components f'jakobsso n et al.] 
l2005iy However, the typical variability of FBQ 0951+2635 
on a timescale of ~ 2 weeks gives us the typical amplitude 
of the deviations S = niB — tua — AmsA, so a correction 
for simultaneity can be estimated from variability studies. 
Throughout most of this paper, we derive Am_BA values 
from single-epoch magnitude differences, whose photomet- 
ric uncertainties are properly enlarged to incorporate the si- 
multaneity error a aim (i.e., the typical amplitude of 8). We 
adopt (Jaim = 0.03 mag, which is consistent with the _R-band 
root-mean-square (rms) variability of FBQ 0951+2635 over 
a 16-day period (see details in subsection 2.1), as well as 
with the i-band variability over such timescale (see subsec- 
tion 3.1). All single-epoch measurements of AmsA include 
the 0.03-mag uncertainty added in quadrature to the pho- 
tometric errors. The flux ratio between components is given 
by A/B = ioO-4x'^'"B-4_ 

2.1 Maidanak follow-up observations 

The Maidanak gravitational lens monitoring programme is 
being conducted by an international collaboration of as- 
tronomers from Russia, Ukraine, Uzbekistan, and other 
countries. The median FWHM (~ 0'.'7) and number of 
clear nights (^ 200 n i ghts per year) at Mt. M aidanak 
(| Artamonov et al] [19871 : lEhgamberdiev et al] [2OO0I ) permit 
to o btain high-resolution images of compact lens systems 
(e.g.. iKoptelova et al.|[2005i ). and here we present and anal- 
yse 7?-band observations of FBQ 0951+2635. These Maid- 
anak homogeneous observations from April 2001 to May 
2006 are an important tool to understand the microlensing 
variability in the double quasar. 

Our Maidanak monitoring consisted of 190 frames (ex- 
posures) in the R band, which were taken with the 1.5m 
AZT-22 Telescope at Mt. Maidanalfl (Uzbekistan) on 37 
different nights (see Table [l]). We used the LN2-cooled 
CCD camera (BroCam) with SITe ST-OOA CCD chip. This 
2000x800 CCD detector has pixels with a physical size of 

^ Information on the Maidanak Observatory is available from the 
Website |http : / / www, astrin . uzsci . net /index, html | 
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Figure 1. Central region of a Maidanak i?-band frame of FBQ 
0951+2635. This frame was taken under very good seeing condi- 
tions (FWHM = 0'.'75) on January 22, 2004. The exposure time 
was 180 s, and we show the logarithm of counts in each pixel. 
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15 /xm, giving angular scales of 0'.'135 pixel" ^ (long- focus 
mode) and 0'.'268 pixel"^ (short-focus mode). The gain and 
readout noise are 1.2 e~ ADU~^ and 5.3 e~, respectively. 

The pre-processing of each frame consists of bias sub- 
traction, overscan trimming, fiat fielding, and cosmic rays 
cleaning. For each observation night, we have two or more 
individual frames obtained under good seeing conditions: 
(FWHM) ~ 1". For example, in Fig.[T]we display the central 
region (I'xl') of an 180-s exposure in subarcsecond seeing 
conditions. This J?-band image (logarithm of counts over the 
FBQ 0951-1-2635 field) includes the SI star (~ 16.6 mag), 
the two components A and B (~ 17.1 mag and ~ 18.3 mag, 
respectively), and the S3 star (~ 19.5 mag). Moreover, the 
average signal-to- noise ratio (aperture radius of 2", i.e., 2x 
FWHM) of an 18 mag star is (SNR) ~ 50. Thus, most of 
the nightly FWHM values are less than the separation of 
the double quasar and there is significant signal associated 
with the faintest quasar component. 

PSF fitting photometry leads to instrumental fluxes of 
components and field stars. These fluxes are calculated using 
the in strumental photometry pipeline by IShalvapin et al.l 
(|2008l ). We also obtain calibrated light curves (in mag) 
of A, B and S3, as well as single-epoch magnitude differ- 
ences ruB — iTiA- The individual values are then combined 
to form nightly means and standard deviations of means 
(i.e., standard errors). The standard errors of niB — niA 
are properly enlarged because we are interested in magni- 
tude differences at the same emission time (see explanation 
above). The fluxes of A, B and S3 on April 15, 2001 are also 
compared to the corre sponding NOT fluxes at close dates 
l|jakobsson et al.]|2005l ). We do not find offsets in the fiuxes 
of A and S3, but the fiux of B is slightly higher (60 mmag) 
than the A^'OT level around Aprfl 15, 2001. This offset is very 
probably due to a small contamination by the lens galaxy 
light (which is absent in the NOT fluxes), so a 60- mmag 
correction is taken into account when obtaining final fluxes 
of B. Each of the four combined curves (m^, ms, mss, and 
Attiba) incorporates 37 data points. These data points and 
other quantities of interest are presented in Table [2] 

Although the standard errors may be good estimators of 
the nightly photometric errors in niA, mB and mss, changes 
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Figure 2. Maidanak _R-band light curves of FBQ 0951+2635. 
Top: A (filled circles), B — 0.7 mag (filled triangles), and S3 — 
1.2 mag (open circles). Bottom: comparison between the A fluxes 
in the first semester of 2004 (filled circles) and the time- and 
magnitude-shifted brightness record of B in that semester (filled 
triangles). See main text for details. 



in the colour coefficient and/or the possible inhomogeneous 
response over the camera area could produce a substantial 
amount of additi onal noise throughout the several years of 
monitoring (e.g.. IShalvapin et al.li2008l '). We check this pos- 
sibility by comparing the average standard error of msa and 
the standard deviation of mss over the whole monitoring 
period. There is a bias factor of about \/3, so we multiply 
the average standard errors of ttia, m_g and msa by \/3 to 
derive typical photometric errors. These final uncertainties 
are 9 mmag (A), 24 mmag (B) and 28 mmag (S3). The top 
panel in Fig.[2]shows the final brightness records of the dou- 
ble quasar and the S3 star. Filled circles, filled triangles and 
open circles trace the behaviours of A, B (shifted by —0.7 
mag) and S3 (shifted by —1.2 mag), respectively. The al- 
most parallel fading by ~ 0.6 mag of the two components 
indicates the presence of long-timescale intrinsic variability 
(extrinsic variations are discussed below). Around day 2750 
(JD— 2451100), one can see an important scatter in the light 



4 V. N. Shalyapin et al. 



Table 2. i?-band photometry of FBQ 0951+2635 from 2001 to 2006. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


10415 


2 




150 


0.94 


0, 


.03 


56.0 


16.884 


0.001 


17.948 


0.015 


1, 


,063 


0.033 


19.490 


0.013 


20304 


4 


1237.797 


180 


0.96 


0, 


.05 


53.0 


16.990 


0.008 


18.056 


0.030 


1, 


,065 


0.048 


19.462 


0.012 


20311 


2 


1 9A4 7Qn 


240 


0.99 


0, 


,19 


53.0 


16.987 


0.003 


18.098 


0.039 


1, 


,111 


0.052 


19.469 


0.004 


20314 


2 


1 9A7 7R9 

J. Zi'i ( . 1 


180 


0.89 


0, 


,10 


53.0 


16.995 


0.005 


18.104 


0.012 


1, 


,108 


0.031 


19.448 


0.035 


30427 


2 




180 


1.06 


0, 


,12 


38.0 


16.995 


0.000 


18.271 


0.009 


1, 


,276 


0.031 


19.458 


0.030 


30521 


2 


J-UOU. UtJU 


180 


0.89 


0, 


,08 


27.0 


17.028 


0.000 


18.210 


0.014 


1, 


,182 


0.033 


19.406 


0.056 


40116 


3 


1 Q9n QQzL 


180 


1.04 


0, 


,06 


64.0 


17.055 


0.011 


18 266 


0.026 


1, 


,212 


0.048 


19.464 


0.035 


40122 


4 


1 Q9R S7Q 


180 


0.91 


0, 


,11 


81.0 


17.062 


0.003 


18 306 


0.007 


1, 


,244 


0.031 


19.455 


0.011 


40209 


5 


1 Qzi4 S4Q 

± y'i'i. OT:y 


180 


0.97 


0, 


,05 


40.0 


17.072 


0.002 


18.296 


0.014 


1, 


,224 


0.033 


19.474 


0.014 


40213 


4 


± yio. yu ( 


210 


1.00 


0, 


,13 


62.5 


17.047 


0.005 


18.304 


0.009 


1, 


,257 


0.033 


19.462 


0.008 


40219 


4 


± yot. yu'i 


180 


0.99 


0, 


,19 


74.5 


17 021 


0.004 


18 346 


0.009 


1, 


,325 


0.031 


19.460 


0.012 


40222 


4 


± yo t .oo 1. 


180 


1.01 


0, 


,17 


44.5 


17 038 


0.002 


18 328 


0.018 


1, 


,289 


0.036 


19.467 


0.026 


40226 


4 


1QR1 S77 
± yu i-.o t i 


180 


1.05 


0, 


,10 


29.0 


17.059 


0.006 


18.365 


0.013 


1, 


,306 


0.035 


19.483 


0.022 


40227 


4 


1 QR9 SRS 
± yui . ouo 


180 


0.91 


0, 


,09 


80.5 


17.036 


0.003 


18.330 


0.009 


1, 


,294 


0.031 


19.439 


0.007 


40301 


4 


1 QR'^ S7R 


180 


0.93 


0, 


,20 


44.0 


17.046 


0.002 


18 333 


0.013 


1, 


,287 


0.033 


19.455 


0.013 




3 


9004 7zlS 


180 


1.13 


0, 


,11 


63.3 


17 092 


0.008 


18 207 


0.017 


1, 


,115 


0.038 


19.483 


0.013 


40413 


4 


901 9 RJ^zl 


180 


0.93 


0, 


,05 


69.0 


17.097 


0.001 


18.295 


0.005 


1, 


,198 


0.030 


19.488 


0.020 


40504 


4 


909Q R77 


180 


0.97 


0, 


,05 


24.5 


17.100 


0.006 


18.327 


0.011 


1, 


,227 


0.034 


19.504 


0.016 


40511 


4 


90'^R R71 


180 


0.99 


0, 


,08 


63.5 


17.068 


0.005 


18 337 


0.009 


1, 


,270 


0.033 


19.444 


0.012 


40520 


3 


.^u'io. uyo 


180 


1.04 


0, 


,04 


60.0 


17 078 


0.004 


18 302 


0.009 


1, 


,225 


0.032 


19.448 


0.007 


41204 


6 




180 


0.90 


0, 


,09 


24.3 


17.208 


0.005 


18.446 


0.010 


1, 


,239 


0.031 


19.480 


0.029 


41205 


g 


^^'io. yoo 


180 


0.89 


0, 


,14 


22.8 


17.217 


0.004 


18.448 


0.012 


1, 


,231 


0.032 


19.479 


0.012 




6 




180 


0.88 


0, 


,06 


62.0 


17.241 


0.003 


18.469 


0.007 


1, 


,228 


0.031 


19.473 


0.008 


50214 


6 


9'^1 S40 


180 


0.97 


0, 


,05 


62.3 


17.226 


0.003 


18.430 


0.012 


1, 


,205 


0.033 


19.475 


0.013 


50306 


6 


2335.815 


180 


0.97 


0, 


,03 


50.0 


17.251 


0.004 


18.481 


0.006 


1, 


,231 


0.031 


19.478 


0.016 


50414 


4 


2374.726 


180 


0.90 


0, 


,04 


47.5 


17.236 


0.006 


18.459 


0.007 


1, 


,224 


0.032 


19.430 


0.006 


51126 


2 


2601.018 


300 


1.09 


0, 


,03 


65.0 


17.339 


0.014 


18.488 


0.015 


1, 


,150 


0.041 


19.489 


0.017 


51128 


2 


2602.974 


300 


0.92 


0, 


,09 


79.0 


17.322 


0.003 


18.569 


0.019 


1, 


,247 


0.034 


19.420 


0.010 


51201 


2 


2606.027 


300 


0.90 


0, 


,10 


61.0 


17.337 


0.002 


18.576 


0.002 


1, 


,239 


0.030 


19.438 


0.003 


60402 


10 


2727.786 


180 


0.83 


0, 


,15 


80.2 


17.372 


0.002 


18.649 


0.005 


1, 


,278 


0.030 


19.377 


0.006 


60411 


10 


2736.793 


180 


1.11 


0, 


,15 


28.4 


17.511 


0.017 


18.457 


0.036 


0, 


,946 


0.060 


19.481 


0.017 


60413 


10 


2738.767 


180 


1.06 


0, 


,08 


29.0 


17.433 


0.006 


18.588 


0.024 


1, 


,155 


0.042 


19.423 


0.024 


60417 


10 


2742.703 


180 


1.05 


0, 


,13 


82.0 


17.428 


0.003 


18.610 


0.007 


1, 


,182 


0.031 


19.401 


0.005 


60509 


10 


2764.678 


180 


1.04 


0, 


,17 


38.6 


17.451 


0.003 


18.637 


0.011 


1, 


,185 


0.032 


19.429 


0.011 


60511 


10 


2766.681 


180 


1.17 


0, 


,10 


38.2 


17.529 


0.019 


18.422 


0.043 


0, 


,893 


0.069 


19.437 


0.014 


60514 


10 


2769.706 


180 


0.97 


0, 


,16 


26.2 


17.483 


0.009 


18.682 


0.012 


1, 


,198 


0.032 


19.453 


0.024 


60517 


10 


2772.664 


180 


1.14 


0, 


,09 


56.6 


17.472 


0.004 


18.528 


0.011 


1, 


,055 


0.033 


19.416 


0.011 



(1) civil date (ymmdd), (2) number of frames, (3) JD— 2451100, (4) exposure time per frame (s), (5) FWHM of the 
seeing disc ("), (6) ellipticity of the seeing disc, (7) SNR of an 18 mag star (within a 2" aperture radius), (8) flux of 
A (mag), (9) error in flux of A (mag), (10) flux of B (mag), (11) error in flux of B (mag), (12) AnigA (™'5'g)i (13) 
error in Amgj\ (mag), (14) flux of S3 (mag), and (15) error in flux of S3 (mag). 



curves of A and B. This scatter was produced when the cam- 
era (BroCam) reached the end of its lifetime. 

In general, sampling is unsuitable when doing analysis 
of short-timescale (intra-year) variability. However, this kind 
of study is possible in 2004 (see Table [l}, i.e., around day 
2000 in Fig. [2] (top panel). In the bottom panel of Fig. [21 
while the filled circles describe the time evolution of the flux 
of A, the filled triangles represent the time- and magnitude- 
shifted light curve of B. The fluxes of B are advanced in 16 
days (time delay; Ijakobsson et al.l |2005| ) and increased by 
1.246 mag (average magnitude difference between the time- 
shifted record of B and the light curve of A; although bins 
with semi-size of 4 days are used, the average magnitude 
difference is within the 1.24-1.25 mag interval for any bin 
semi-size below 6 days). We show that there is good agree- 
ment between both brightness records, i.e., the two records 
are consistent with each other in the overlap periods. For 
this reason, we can state that short-timescale microlensing 



is absent or elusive. The time-delay-corrected flux ratio [R- 
band data from day 1900 to day 2050) isA/B = 3.15 ± 0.05 
(la). 

We also obtain the Maidanak and NO T J?-band typical 
variability of FBQ 0951-1-2635 at time lags Ai 100 days, 
with special emphasis on At = 16 days, i.e., the time delay 
between components. The typical magnitude variation on 
this timescale is directly related to the simultaneity error 
Cairn that we introduce in the second paragraph of section 
2. The rms variability {(As)'^)^/^ at lag At is given by 

{{MY' = {(1/7V) ^[(m, - m^f - al - (1) 

where the sum only includes the (i,j) pairs verifying that 
t] — ti ^ lS.t (the number of such pairs is A*'). To derive 
the rms values from Eq. (1), we consider the magnitude 
variations in both components A and B (Maidanak fluctu- 
ations from day 2700 to day 2800 are doubtful, and thus. 
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this Maidanak observation period is not taken into account 
in the analysis; see above). Fig. |3] shows the structure func- 
tions obtained from the Maidanak (filled circles) and NOT 
(filled squares) light curves, using independent 4-day bins. 
Two vertical dashed lines define the bin of interest (around 
the lag of 16 days). At lags At < 60 days, the quasar was 
somewhat less active over days 50—1000 (NOT data). The 
horizontal dashed line corresponds the 0.03 mag level, which 
roughly coincides with the Maidanak rms variability at the 
lag centered on 16 days. These results explain why we take 
Csim = 0.03 mag, at least at red wavelengths. We remark 
that the bottom panel of Fig. [2] shows a "flare" of ~ 0.1 mag 
over ~ 10 days, between day 1990 and day 2000, and this 
seems to question our variability analysis in Fig. [3] How- 
ever, taking the error bars into account, the true variation 
in the flux of B could be of about 30 mmag. This is in good 
agreement with the adopted variability over a 16-day period. 
Moreover, the ~ 0.1-mag jump (considering central values) 
is only deflned by two consecutive fluxes of the B compo- 
nent. In other words, it is a very poorly sampled variation 
in the faintest and noisest component, which is most likely 
caused by observational noise. 

The 37 single-epoch measurements of AmsA in Table [2] 
are grouped in 7 different time intervals (see the seven pe- 
riods in Table [T] and the corresponding days in Table |2]). 
We compute the weighted average and its uncertainty for 
each interval. These averages are then translated into 7 R- 
band flux ratio measurements with an 1—3% accuracy. The 
new flux ratio values are depicted in Fig. |4] (filled circles) . 
For example, the flux ratio around day 2000 (3.17 ± 0.03; 
la) is consistent with the result in two paragraphs above, 
which takes the time delay correction into account. Around 
day 2750, we obtain two estimates of A/B due to the fact 
that the photometric data show a large scatter (see above). 
Firstly, we use all available nights (last open circle in Fig.|4]). 
Secondly, only high-quality nights, i.e., FWHM < I'.'l and 
SNR > 50, are considered (filled circle above the last open 
circle in Fig.|4]). This second estimation seems more reliable, 
so the first one is assumed to be biased. We also compare 
the A'^OTflux ratio at an intermediate epoch and the corre- 
sponding Maidanak ratio from the former ST-7 CCD camera 
(which was operating in that epoch). The Maidanak/ST-7 
1(7 measurement, 2.58 ± 0.05 (first open circle at day 484), 
is in excellent agreement with the A'^OT gradient and limits 
(dashed lines; see the first paragraph in subsection 2.3). 
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Figure 3. i?-band rms variability at different time lags. The 
Maidanak (filled circles) and A'^OT (filled squares) trends are anal- 
ysed in a separate way. We use 4-day independent bins centered 
on lags of 4, 8, 12, 16, 20,... days. The bin of interest (around the 
lag of 16 days) is highlighted from two vertical dashed lines, and 
the horizontal dashed line represents the 0.03 mag level. 



r-band flux of B is not corrected by a possible contami- 
nation by the lens galaxy light (the r-band correction will 
be significantly smaller than the 60-mmag offset in the R 
band, since the r Sloan filter does not transmit light at 
wavelengths of 700-900 nm), and the r-band magnitude 
difference is taken as the magnitude difference in the R 
band (this is a reasonable assumption because the r — R 
colours of both components are similar). We note the 
agreement between the SDSS flux ratio at day 2250 (flrst 
flUed triangle in Fig. |4]), and the Maidanak/BroCam value 
of A/B around day 2300 (filled circle next to the flrst flUed 
triangle in Fig. |4} . 

Maidanak observations in 2007 were done with the 
Spectral Instruments (SI) 600 Series CCD camera. It is a 
4096x4096, 15-fim pixel CCD. The rectangular pixels have 
angular scales of 0'.'303 (along the horizontal axis) and 0'.'258 
(along the vertical axis). The camera gain is 1.45 e~ ADU~^ 
and the readout noise is 4.7 e~. These Maidanak/SI 600 Se- 
ries frames belong to a private archive that is managed by 
several teams. To test the time evolution of A/B in recent 
years, we use the three _R-band frames that were taken on 
April 13, 2007. After applying a 60-mmag correction to the 
fluxes of B, the flux ratio at day 3104 is A/B = 3.17 ± 0.10 
(last filled triangle in Fig. H]). 



2.2 Complementary frames 

A few complementary frames are used to check 
the reliability of the A/B values from the Maid- 
anak/Bro Cam homogeneous monitori ng. The SDSS 
archiv^fl l|Adelman-McCarthv et al]|2008l ) contains a frame 
of FBQ 0951-1-2635 in the r band. This frame was taken 
under very good seeing conditions (FWHM = 0'.'83) on 
December 11, 2004. From PSF fitting photometry, we infer 
m_B — mA = 1.22 ± 0.01 mag (Icr), which translates into 
AniBA = 1.22 ± 0.03 mag and A/B = 3.08 ± 0.08. The 

2 See the SDSS Data Release 6 (DR6), which in- 
corporates the DR6 Catalogue Archive Server site at 
|http://cas.sdss.org/ astrodr6 /en/ 1 



2.3 Global perspective on flux ratio evolution 

IParaflcz et al] l|2006h suggested two possible 7?-band mi- 
crolensing gradients over day s 50—1000 (see mid dle and 
bottom panels in Fig. 2 of iParaficz" et al.1 l2006l ). How- 
ever, the second solution (0.077 ± 0.00 7 mmag day~^) 
seems more consistent with the analysis bv lJakobsson et aU 
(l2005l). so we only con sider this early microlensing slope. 
IJakobsson et al.l (2005h found a flux ratio A/B ^ 2.65 ± 
0.02 (Am_BA ~ 1.06 ± 0.01 mag) when they focussed on the 
last part of the NOT brightness records, i.e., around day 
800 (see the filled square in Fig. [4l). This is fully consistent 
with the second linear fit bv lParaflcz et all (|2006l ) (see bot- 
tom panel in Fig. 2 of that work), which indicates AmsA 
= 1.06 at day 800. Alternatively, the extrinsic variability 
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Figure 4. Time evolution of the optical continuum (red wave- 
lengths) flux ratio of FBQ 0951+2635. We show a previous mea- 
surement from NO T i?-band frames (filled square at day 800) , as 
well as new Maidanak/BroCam _R-band estimates (filled circles) 
and two additional results from SDSS and Maidanak/SI 600 Se- 
ries frames (filled triangles). We also show the lower limit of the 
Mgll flux ratio (base of the vertical arrow). The open circles and 
dashed lines are explained in the main text. 



could be due to two consecutive short-term microlensing 
events before day 700. While the first event was fitted to 
a Gaussian, it is also apparent the possible existence of a 
second event (see data point distribution between days 400 
and 600). From day 700 to day 1000, the observations are 
consistent with constant microlensing because almost all er- 
ror bars cross the zero level. If the two rapid events with 
similar amplitude are true features, no gradient is required 
to account for the observations. The two alternatives and 
their extrapolations to more recent epochs are depicted in 
Fig. |4] The inclined dashed lines represent the gradient of 
0.077 ± 0.007 mmag day~^ (absence of short-term events). 
The horizontal dashed lines are the lower and upper limits 
associated with the second scenario: absence of a gradient, 
but presence of short-timescale fluctuations with amplitude 
of about 60 mmag (flaring behaviour of B). 

As we can see in Fig. |3J the Maidanak measurements 
from day 2600 to day 3100, i.e., the last two filled circles and 
the last filled triangle, are consistent with the first microlens- 
ing scenario. However, the most recent behaviour oi A/B 
does not imply the presence of a long-term gradient lasting 
~ 10 years. The whole set of results within the ~ 3000-day 
monitoring period suggests the existence of a long-timescale 
microlensing fluctuation, which had a bump between day 
1200 and day 2200, and has remained almost flat during 
the last years. The bump consists of a signiflcant increase 
(in A/B) of about 15% followed by a shallow decrease. As 
the flux ratio is measured accurately, the signal-to-noise ra- 
tio for the prominent increase is ~ 3—5. When a source 
crosses a microlensing magniflcation pattern, the associated 
light curve may have a complex structure including different 
gradients at different epochs (e.g.. IWambs ganss 199^). We 
also point out that our new flux ratio estimates clearly dis- 
agree with the second scenario (see horizontal dashed lines 
in Fig. |4ll, which is based on the exclusive production of 
short-term (~ months) fluctuations. 



3 FOLLOW-UP IN THE i SLOAN PASSBAND 
3.1 LQLM campaign 

We included FBQ 0951-1-2635 as a key target in our 
Liverpool Quasar Lens Monitoring (LQLM) programme 
iGoicoechea et all l2008l : IShalvapin et all 120081 ) . The opti- 
cal frames were taken with the i-R2Q between February 6, 
2007, and May 31, 2007. Each observation night consisted 
of five exposures of 100 sec in the i band, using a dither 
cross pattern (see observations summary in Table [TJ . The 
ii?r pre-processing pipeline included bias subtraction, trim- 
ming of the overscan regions, and flat fielding. These pre- 
processed frames are publicly available on the Lens Image 
Archive of the German Astrophysical Virtual Observatorjy. 
We initially remove all frames that either are characterised 
by an anomalous image formation or have large seeing discs 
(FWHM > 2" in the frame headers). Later we carry out 
cosmic rays cleaning and defringing. In the last step of the 
whole pre-processing procedure, all frames in each night (^ 
5) are combined, i.e., they are aligned and then averaged. 

We create a stacked frame (consisting of the combina- 
tion of the best exposures) to better detect the lens galaxy 
and to subtract its light in the combined frames. This repre- 
sents a total exposure time of 4.4 hours. The stacked image 
is characterised by FWHM = 1'.'17 and a very high signal-to- 
noise ratio (SNR = 413), where FWHM and SNR (aperture 
radius of 2x FWHM) are measured on the SI and S3 stars, 
respectively. In the i band, S3 has a magnitude similar to 
that of the faintest quasar component B. Unfortunately, nei- 
ther our best combined frames in terms of FWHM and SNR, 
nor our deep stacked frame lead to detection of the lens 
galaxy. Moreover, we obtain meaningless results by using 
constraints on the rela tive position and brigh tness profile of 
the galaxy (taken from ljakobsson et al.l2005l ) . The 2D signal 
from the very faint galaxy seems to be strongly affected by 
the photon noise of the quasar components and background, 
so we cannot resolve that signal in our frames. Th us, our in- 
strumental photometry pipeline jShalvapin et al. I2OO8) only 
fits the instrumental fluxes of both components, which are 
described by stellar-like objects (PSF fitting photometry, us- 
ing SI as PSF star). Once the instrumental photometry is 
done, we obtain calibrated and corrected brightness records 
of A, B and field stars. The transformation (calibration- 
correction) pipeline incorporates zero-point, colour and lin- 
ear inhomogeneity terms, and the fin al magnitudes are give n 
in the SDSS photometric system (jShalvapin et al.l |2008| ). 
This transformation software is exclusively applied to the 
combined frames with FWHM < l'.'5 and SNR > 50, i.e., on 
22 out of 52 observation nights (see Tabled. 

The top panel in Fig. [S] displays the LQLM light curves 
of the lensed quasar (A and B components) and the S3 star. 
Filled circles, flUed triangles and open circles trace the be- 
haviours of A, B (shifted by —0.9 mag) and S3 (shifted by 
—0.6 mag), respectively. There must be some diffuse con- 
tamination by the lens galaxy light, so the B fluxes in Table 
3 and Fig. [S] are greater than the true ones. From the stan- 
dard deviation of mss over the whole monitoring period, we 



^ Information on the Liverpool Robotic Telescope is available 

from the Website http://telescope.livjm.ac.uk/j 

^ See the Web site h ttp://vo.uni-hd.de/ lensdemo / view /q/form | 
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Table 3. j-band photometry of FBQ 0951+2635 in 2007. 
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(1) civil date (ymmdd), (2) number of 100-s exposures, (3) JD-2451100, (4) FWHM of the seeing disc {"), (5) 
ellipticity of the seeing disc, (6) SNR of the S3 star (within a 2x FWHM aperture radius), (7) flux of A (mag), (8) 
flux of B (mag), and (9) flux of S3 (mag). 

infer a typical error in the S3 fluxes of 17 mmag (see open 
circles and associated error bars). Hence, the typical uncer- 
tainty in the B fluxes should be larger than 17 mmag, since 
B is as bright as S3 (~ 18.5—18.6 mag), but it is located 
within a crowded region. Unfortunately, even using frames 
with (FWHM) ~ 1'.'2 and (SNR) ~ 80, we cannot achieve 
1—2% level photometry for B (some deviations between ad- 
jacent nights are relatively large). 

Four pairs of adjacent fluxes of A show significant scat- 
ters around the mean values (in the top panel in Fig. [S] 
we draw four lines joining the members of each pair). For 
this reason, they are grouped (by computing four mean 
fluxes at the corresponding mean epochs) to obtain a mas- 
ter light curve, i.e., an accurate (smooth) trend that reli- 
ably describes the underlaying short-timescale variability. 
This master curve is depicted in the bottom panel of Fig. [S] 
(filled circles). Before the gap, the quasar has a small level 
of activity. However, after the gap, there is a 60-mmag gra- 
dient lasting about 30 days. Thus, the _R-band and i-band 
records exhibit similar variability levels on a timescale of ~ 
2 weeks (see Fig. [3]). To gain perspective on the nature of 
the observed variability in the i band (A component), the 
B light curve is shifted in ti me (using the 16-day time de- 
lay by Jakob sson et al.ll2005l ). and then is compared to the 
master light curve of A. We derive an average magnitude 
difference of 1.094 mag between the time-shifted record of B 
and the master curve of A (using bins with semi-size of 2-4 
days). The time- and magnitude-shifted light curve of B is 
also included in the bottom panel of Fig. [5] (filled triangles) . 
The typical error in the S3 fluxes is used as a lower limit 
for the uncertainty in the photometric measurements of B 
(see above). We flnd a clear agreement between the master 
curve of A and the adjacent fluxes of B, when these last 
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Figure 5. LQLM i-band light curves of FBQ 0951-1-2635. Top: A 
(filled circles), B — 0.9 mag (filled triangles), and S3 — 0.6 mag 
(open circles). Bottom: comparison between the A master record 
(filled circles) and the time- and magnitude-shifted light curve of 
B (filled triangles). See main text for details. 
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fluxes are properly shifted in time and magnitude. Thus, 
there is no evidence of short-timescale extrinsic variabihty 
(due to microlensing or another phenomenon) in our i-band 
observations. 

The new LQLM hght curves of FBQ 0951+2635 aUow 
us to measure an accurate flux ratio (corrected by the time 
delay) in the i band: ^/S* = 2.74 ± 0.02 {la). We do not 
measure the i-band flux ratio A/B around day 3100 but the 
contaminated ratio A/B,, where B, = B + G and G is a 
contribution due to the lensing galaxy. 



3.2 SDSS frame and recent evolution of A/B 

The high-quality SDSS^ i-band frame of FBQ 0951-^2635 at 
day 2250 (FWHM = 0'.'75 and SNR = 99) leads to a flux 
ratio A/B-j, — 2.8 ± 0.1 (If), in good agreement with the 
LQLM estimation. We note that B, = B + G (see above), 
and the uncertainty in the SDSS flux ratio incorporates both 
photometric and simultaneity errors. This last error is asso- 
ciated with the use of fluxes at the same time of observa- 
tion, i.e., without time delay correction (see above). There is 
no evidence of an appreciable evolution in A/B, over days 
2250—3100. This i-band result agrees with the Maidanak- 
SDSS recent trend of A/B in the R band (see filled circles 
and triangles from day 2250 to day 3100 in Fig. |4]). Com- 
paring the flux ratios in both optical fllters, we may derive 
that the lens galaxy light in the i band produces a plausi- 
ble contamination of the B component in such fllter (~ 130 
mmag) . 



4 CONCLUSIONS 



A previous analysis (jParaficz et al.l [20061 ) showed the exis- 
tence of extrinsic variability of the observer-frame optical 
continuum in the gravitationally lensed double quasar FBQ 
0951-1-2635. The two quasar components (A and B) cross two 
different regions of the lensing galaxy, so distributions of col- 
lapsed objects could affect one (or both) of the light curves 
and then produce extrinsic variations. This microlensing hy- 
pothesis is more plausible for the B comp onent because it 
cross es the central reg i on of the galaxy (see lJakobsson et al.l 
bOOSl ). iParaflcz et al.l (|2006l ) reported on two possible mi- 
crolensing variabilities at red wavelengths in the 1999—2001 
period: short-timescale events (B flares having a duration of 
months) or a gradient lasting a few years. Both alternatives 
can account for the 1999—2001 _R-band observations, and 
this paper sheds light on the dominant kind of microlensing 
fluctuations occurring in FBQ 0951+2635. 

We analyse new Maidanak _R-band images taken in the 
2001—2006 period (a 6-year homogeneous monitoring), and 
a few complementary frames in the rR bands taken in 2004 
and 2007 (SDSS and Maidanak/SI 600 Series archives). The 
Maidanak-SDSS flux ratios (A/B) in the red part of the 
optical continuum are inconsistent with the absence of long- 
timescale (~ years) gradients and the continuous production 
of B flares lasting a few months (short-timescale microlens- 
ing events). If this last scenario were true, all data points 
in Fig. |4] would be distributed between the two horizontal 
dashed lines. The whole set of flux ratios (from 1999 to 2007) 
favour the existence of a long-timescale microlensing fluctu- 
ation, so long-timescale gradients seem to be the dominant 



microlensing variations. While A/B shows a bump in the 
2003—2004 period, it is almost constant from late 2004 to 
the middle of 2007. This last quasi-stationary behaviour is 
supported by additional data in the i band (from new LQLM 
ima ges taken in 2007 and the SDS S archive frame in th at fil- 
ter). [c^nuliinaeEaLl (|2005f ) and lParaficz et al.l (120061 ) also 
found a long-timescale microlensing variation in the dou- 
bly imaged quasar SBS 1520+530. Although the existence 
of rapid fiares in the early years of monitoring (1999—2001) 
cannot be ruled out, these hypothetical B fiares are not pro- 
duced in a continuous way. Short-timescale microlensing is 
not detected in the Maidanak variability study over the first 
semester of 2004. The LQLM data in the i band also indicate 
the absence of short-timescale events in 2007. Apart from the 
1999-2007 optical d ata, FBQ 0951+2635 ha s recently been 
observed in X-rays l|Dai fc KochanekllioogI '). However, the 
X-ray flux ratio has a large uncertainty and it is not useful 
for comparison with our estimates in the optical continuum. 

If microlensing in the B component is the physical ori- 
gin of the long-timescale fluctuation, this fluctuation rep- 
resents a progressive demagnification of B followed by a 
quasi-stationary evolution in more recent epochs. The opti- 
cal continuum fiux r atio at red wavelength s is always below 
the Mgll flux ratio l|jakobsson et ah I l2005l . see also Fig.©, 
which is usually assumed to be weakly affected by microlens- 
ing. Thus, the recent microlensing magnification of B would 
be relatively weak, but still appreciable. The microlensing 
peak (maximum magnification of B) would have taken place 
before the discovery of the lens system by ISchechter et al.l 
(1998). Strictly speaking, the continuum fiux ratio refers to 
the flux ratio in the R band, i.e., a spectral region containing 
both the red continuum and the Mgll emission line. How- 
ever, this line only contri butes a few percent to the broad 
band fluxes (see Fig. 4 of ISchechter et al]|l998l ). It can be 
shown that the line-corrected ratio is very similar to the ra- 
tio fr om total fluxes in the R band (see also I Jakobsson et al.l 
l2005h . We also note that A/B is not corrected by any ex- 
tinction factor, so this fiux ratio might be a biased estimator 
of the lens (m acrolens + microlens) magnification ratio. For 
example, Fal co et al.l ((l999) suggested the relative extinc- 
tion of A by dust in the lensing galaxy. This dusty scenario 
leads to a true lens magnification ratio greater than A/B. 

We point out that new flux ratios at a large 
collection of wavelengths should be key tools to dis- 
cuss dust ex tinction and microlensing in the lensing 
galax y ("e.g.. [ Wambganss fc Paczvnskil 199ll: iFalco et al.l 



19991: iGoicoechea et al.l l2005al lbl: lEh'asdottir et al.1 \2Q(M . 



Yonehara et al.ll2008l ). Moreover. FBQ 0951+2635 must be 



imaged in the rR bands during the next years (several times 
per year, with a 2-week separation between consecutive ob- 
servations). This modest monitoring programme will lead 
to draw the future evolution of A/B, as well as to obtain 
relevant information on the structure of both the source 
quasar and the intervening galaxy ( e.g., I WambsganssI 1 1990l : 
lKochanek|[2004l : iMorgan" et al. I l2007l ). 
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